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Scale Determination Using the Static Potential with Two Dynamical 
Quark Flavors 
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Indiana University, Bloomington, IN 47405, USA 

We study the static potential using gauge configurations that include the effects of two flavors of dynamical 
Kogut-Susskind quarks. The configurations, generated by the MILC collaboration, and used to study the spectrum 
and heavy-light decay constants, cover a range 5.3 < < 5.6. There are at least four quark masses for each 
coupling studied. Determination of ro from the potential can be used to set a scale. This alternative scale is 
useful to study systematic errors on the spectrum and decay constants. 



1. INTRODUCTION 

To get an idea of systematic errors introduced 
in determining the physical scale in lattice QCD 
simulations using nip, it is important to use an al- 
ternative quantity to set the physical scale. One 
quantity, tq, proposed by Sommer |1[ has recently 
been studied in several papers (2|-|||7 Here, we ex- 
amine the static potential for a range of couplings 
5.3 < 6/(7^ < 5.6 with two flavors of dynamical 
staggered quarks. We compute the Sommer scale 
for each dynamical quark mass and extrapolate to 
the chiral limit at each coupling. Finally, using 
the scale determined by ro, we take the contin- 
uum limit and replot the Edinburgh curve, find- 
ing very little difference as compared with setting 
the scale from mp. 

The configurations were generated by the 
MILC collaboration |^ using the staggered quark 
action with two dynamical flavors and the Wilson 
gauge action. They were stored every 10 units of 
molecular dynamics time. These configurations 
have been used elsewhere for light hadron spec- 
trum calculations For each coupling there 
are at least four quark masses. (For details, see 
Ref. §.) 

2. STATIC POTENTIAL 

To calculate the static potential, Wilson loops 
W{R, T) were measured for 5 and 10 smearing 
steps using the APE smearing method as ex- 
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Figure 1. Hint of string breaking: flattening of 
the last few points. (3 = 5.3, amq = 0.075. 



plained in The effective potentials were then 
calculated using 



Vt{R) = In 



W{R,T) 



(1) 
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W{R,T+l) 

After blocking the results into 30 time unit blocks, 
errors are obtained using single elimination jack- 
knife. V{R) is expected to flatten out at large R 
for dynamical quarks due to string breaking. At 
our strongest coupling, 5.3, and smallest quark 
masses we have searched for evidence of string 
breaking. For the two lightest masses, we did not 
find a signal at large enough values of i?, but we 
do find a hint of string breaking for avriq — 0.075. 
In Fig. 1, the potential seems to flatten for R 
greater than five to seven lattice spacings. This 
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Figure 2. The graphs show ro/a vs. arUq for 
different /3 values, along with an exponential fit 
defined in Eq. (4). 

corresponds to a distance of 1.3-f .8 fm. Ref. |^ 
reports finding signs of string breaking at 0.8-1. f 
fm, however their mTr/wp is lower. This is con- 
sistent: the string is expected to break at smaller 
distances for lighter quarks. 
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Figure 3. Continuum extrapolation for m^/mp = 
0.1753. Continuum extrapolation was made at 
various m^/mp values to plot the Edinburgh 
curve shown in Fig. 4 

3. FIT PARAMETERS 

We used the standard ansatz for the potential, 
V{R) = Vo + (TR-e/R- f{GL-l/R). (2) 

This is a simple combination of a Coulomb term 
at short distances and a linear increase at large 
distances. The last term accounts for lattice ar- 
tifacts. For more details on this ansatz and its 
limitations see Refs. § and 0. 

Using the fit parameters a and e, the Sommer 
scale ro is calculated in the usual way, 

-»=V^ (3) 

with c = 1.65. This corresponds to a physical 
length of tq ^ 0.5 fm. 

We fit the potential to the ansatz (2) for dif- 
ferent ranges of R and different values of T. The 
" good" fits yielded consistent values, except for (3 
— 5.35, amg=0.05. We have excluded this point 
from our fits for the chiral extrapolation. The 
criteria for best fits were the confidence of the 
fit and the degrees of freedom. Additionally, the 
range of R had to include the value of tq, which 
was an important constraint for the coarse lat- 
tices. In Fig. 2, we show the best estimates of tq 
along with the exponential fit explained below. 
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4. SCALE DETERMINATION 



the scale determined here in Fig. 4. 



For dynamical quarks, a different value of tq 
and a is obtained for each quark mass. The con- 
sideration that the physical scale should approach 
the quenched result as the quark mass tends to 
infinity and the graphs of ro vs aniq suggest an 
exponential fit 



ra = A + B\ 



1) 



(4) 



with three parameters A, B and C. 

We obtained good fits to this form for all /? val- 
ues except 5.415, where we excluded the heaviest 
two points from the fit to get an acceptable con- 
fidence level. These fits were used for the chiral 
extrapolation to determine ro for each niq. 
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Figure 4. Continuum Edinburgh curves for dy- 
namical fermions using the scale as determined 
by TJip and ro. 

5. EDINBURGH PLOT 

A most surprising result of Lattice '98 was 
that the Edinburgh plot for dynamical quarks was 
shifted up compared to the quenched simulations, 
leading to a larger deviation from the real world 
value in the continuum limit We were curi- 
ous to know if the Edinburgh curve changed when 
setting the scale from Tq. Fig. 3 shows mjv/TOp 
as a function of a/r^ for the real world mT^jmp 
value .0 The Edinburgh curve is replotted using 

■^These fits were good up to m-n:/mp = 0.45 but for higher 
ratios the /3 = 5.415 point departed from the curve. There 



6. CONCLUSIONS 

We have calculated the heavy quark potential 
for a range of couplings for dynamical staggered 
quarks. Some flattening of the static potential 
for large R at strong coupling is seen which hints 
at string breaking, but does not suffice to make 
a conclusive case. We determine rg, extrapolate 
it to the chiral limit using an exponential form 
and examine the hadron mass data using the new 
scale. We do not find significant deviations from 
the previous continuum extrapolations using the 
p mass to set the scale. 
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Origin-2G00, a Paragon and the CANDYCANE 
Linux cluster at Indiana University. 

REFERENCES 



R. Sommer, |hep-lat/9310022|. 

R.G Edwards et al, |hep-lat/9711003| v2. 
U.M. Heller et al, Phys. Lett. B335, (1994) 
71. 



C.R. AUton et al, |hep-lat/9808016| , and ref- 
erences therein. 

C. Bernard et al, Nucl. Phys. B(Proc. Suppl.) 
60A (1998) 297. 



9. 

10. 



C. Bernard et al, |hep-lat/9810035 . 

S. Tamhankar and S. Gottlieb, in preparation. 

M. Albanese et al, Phys. Lett. 192B (1987) 

163. 

C. DeTar et al, these proceedings. 

C. Michael, Phys. Lett. 283B (1992) 103. 



seems to be a problem with that couphng for higher masses 
that we do not completely understand yet. Hence we have 
excluded that point from the fits for continuum extrapo- 
lation for nijr/'m.p > 0.45. This exclusion yielded better 
fits without much shift in the continuum value. 



